AnEnergyConserving Local DiscontinuousGalerkin Method for a Nonlinear Variational Wave Equation

We design and numerically validate a local discontinuous Galerkin (LDG) method to compute solutions to the initial value problem for a nonlinear variational wave equation originally proposed tomodel liquid crystals. For the semi-discrete LDG formulation with a class of alternating numerical fluxes, the energy conserving property is verified. A dissipative scheme is also introduced by locally imposing some numerical “damping” in the scheme so to suppress some numerical oscillations near solution singularities. Extensive numerical experiments are presented to validate and illustrate the effectiveness of the numerical methods. Optimal convergence in L2 is numerically obtained when using alternating numerical fluxes. When using the central numerical flux, only sub-optimal convergence is observed for polynomials of odd degree. Numerical simulations with long time integration indicate that the energy conserving property is crucial for accurately capturing the underlying wave shapes. AMS subject classifications: 65M60, 65M12, 35Q35